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SUMMARY 

A diagnostics facility for MPD thruster plume measurements has been built 
and is currently undergoing testing. The facility includes electrostatic 
probes for electron temperature and density measurements. Hall probes for mag- 
netic field and current distribution mapping, and an imaging system to estab- 
lish the global distribution of plasma species. Preliminary results for MPD 
thrusters operated at power levels between 30 and 60 kW with solenoidal applied 
magnetic fields show that the electron density decreases exponentially from 
1x1 0^0 m - 3 to 2 x 10^8 m~3 over the first 30 cm of the expansion while the elec- 
tron temperature distribution is relatively uniform, decreasing from ~2.5 to 
1.5 eV over the same distance. The radiant intensity of the Aril 4879 A line 
emission also decays exponentially. Current distribution measurements indicate 
that a significant fraction of the discharge current is blown into the plume 
region, and that its distribution depends on the magnitudes of both the dis- 
charge current and the applied magnetic field. 


INTRODUCTION 

The potential for high specific impulse and system simplicity make MPD 
thrusters attractive for propulsion applications on planetary, orbit transfer, 
and maneuvering missions (refs. 1 and 2). However, demonstrated thruster effi- 
ciencies are still too low to make MPD thrusters practical for these applica- 
tions. The goal of this research is to establish the performance limitations 
of these devices by examining the physical mechanisms of thruster energy 
deposition. 

The physical processes which govern acceleration and internal energy depo- 
sition via ionization and propellant heating are distributed throughout the 
thruster and near field plume. Interelectrode electron densities measured in 
quasi-steady thruster (1 to 2 ms pulse) by Turchi (ref. 3) and Tahara (ref. 4) 
indicate that the primary ionization zone for argon propellant is near the 
thruster backplate. However, Bruckner (ref. 5) and Boyle (ref. 6 ), using opti- 
cal and double electrostatic probe techniques, measured the exhaust velocity of 
self-field thrusters and found that significant acceleration took place in the 
first 15 cm outside the thruster. Maisenhalder (ref. 7) used rapidly rotating 
single electrostatic probes to establish flow streamlines, electron tempera- 
tures, and densities in the plume of a steady-state, self-field MPD thruster. 
Fradkin (ref. 8 ), using a Rogowski coil, showed that as much as 40 percent of 
the discharge current could flow more than 90 cm downstream of an applied field 
MPD thruster. These measurements indicate a substantial amount of energy depo- 
sition in the plume region, but are limited to high power quasi-steady or low 
power steady-state operation and to a narrow range of thruster operating 
conditions. 
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The goal of this work is to establish the plasma characteristics and domi- 
nant physical processes present in the exhaust plume of both self and applied 
field steady-state MPD thrusters operated at power levels between 20 and 
250 kW. Following an overview of the experimental facility, the plume diagnos- 
tic techniques and their implementation are presented. Preliminary measure- 
ments of the electron density and temperature, ionized argon line intensity and 
current distributions are shown and their implications discussed. 


EXPERIMENTAL FACILITY AND THRUSTER CONFIGURATIONS 

The MPD thruster test facility consisted of a 3 m diameter, 3 m long spool 
piece attached to a 21 m long, 7.6 m diameter vacuum tank via a 3 m diameter 
gate valve. Thrusters were mounted on a thrust stand within the spool piece, 
and the gate valve was opened to test the thrusters. The thruster support 
assembly is shown schematically in figure 1. The tank pressure was maintained 
at pressures below 0.07 Pa <5x1 0 -4 torr) during testing at argon flow rates up 
to 0.16 g/s. The thrusters were powered by four 66 kW welding supplies com- 
bined in a series/parallel circuit to give 88 V at 3000 A. A solenoidal mag- 
netic field was applied using a 46 turn, 14 cm inner diameter, 15.2 cm long 
coil surrounding the thruster anode. A single welding supply was used to drive 
up to 1500 A through this magnet, generating a 0.2 tesla magnetic field at the 
cathode tip. Further details on the facility may be found in references 9 
to 11. 

The thruster configurations used for these experiments, shown in figure 2, 
used a cylindrical water cooled copper anode and a 2 percent thoriated tungsten 
cathode. The thruster designations (D, E, and F) correspond to those used in 
reference 10. Configuration D had a cylindrical anode which ended in a short, 
41° flare. The continuously flared anode geometry used for geometries E and F 
was 5.4 cm long with a 10.5° half angle. For all thrusters the cathode diame- 
ter was 1.27 cm; its length was 2.2 cm for geometries D and E and 3.8 cm for 
geometry F. The primary chamber dimensions are listed in table I. Argon pro- 
pellant was injected through the backplate. 


DIAGNOSTICS 

Electrostatic and Hall probes were used with a positioning system to study 
the electron density/temperature and current distributions and a charge injec- 
tion device (CID) camera was used with narrow bandpass interference filters to 
establish global plasma species distributions. In this section detailed 
descriptions of these diagnostics are presented. 


Probe System 

The principal problem with intrusive diagnostics of high power steady- 
state plasma thrusters is the extremely hostile environment in which the meas- 
urements are made. Total enthalpies in the exhaust plume can reach up to 
5x10^ J/kg. For this research the probes were swept through the plume in 
~50 ms to avoid the large sizes associated with active cooling. The drive 
system, shown in figure 3, consisted of a stepper motor driven linear actuator 
used to move the probes axially on which a dc motor was mounted to rotate the 
probes through the plume. The probe arm length was such that the probes passed 
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directly through the thruster axis and a variable potentiometer was used to 
detect the angular position of the probe arm during the plume traverse. Data 
were acquired at 1 kHz during each pass through the plume. The entire motion 
and data acquisition sequence was controlled using a personal computer. 

Electrostatic probes . - The rapid probe motion makes use of standard 
single and double probe techniques (ref. 12) difficult. For this reason the 
simpler triple probe technique developed by Chen (ref. 13), and recently 
implemented at the University of Stuttgart (ref. 14), was used for these exper- 
iments. This technique eliminates the requirement for a voltage ramp and per- 
mits rapid evaluation of the electron temperature and density. Comparisons of 
results from the triple probe method with double probe and nonintrusive inter- 
ferometric techniques have shown good agreement under conditions similar to 
those expected in MPD thruster plumes (refs. 15 and 16). As shown in figure 4, 
the triple probe technique relies on applying a constant bias voltage between 
two probes while a third is allowed to electrically float nearby. The electron 
temperature, T e , and density, N e , are related to the floating voltage, V^, and 
the current, 1 ^ 3 , by (ref. 13) 
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where the voltages are defined in figure 4 and Ap is the individual probe 
surface area. These relations are based on the assumptions of a Maxwellian 
electron distribution function, a collisionless probe sheath, and negligible 
inter-probe interaction. The latter requires that the probe separation be much 
larger than the Debye length. The probes are aligned with the thruster axis to 
minimize the impact of flow angle on the electron density measurements. This 
effect goes to zero along the axis (where the probes are in line with the 
flow), and the axial variation in the shape of the electron density distribu- 
tion implies this effect introduces an error of less than 20 percent at the 
largest axial/radial position examined (worst case). The impact of the mag- 
netic field on the probe measurements should be small for the conditions 
studied (ref. 13). 


The triple electrostatic probe consisted of three 0.25 mm diameter, 1 cm 
long, tungsten wires each separated by 1 mm. Probe end effects (ref. 17) 
should be small for this geometry (length to diameter ratio, "50) . Alumina 
tubes, which were 1.4 mm in diameter and 3.8 cm long, were used as the insula- 
tors. A constant bias (V^) of 17 V was applied using a 30 V, 2 A power sup- 
ply. Testing revealed that varying between 6 and 18 V had little impact 

on the measurements. All signals were sent through isolation amplifiers (3 dB 
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point at 10 kHz) before going to the data acquisition system. Repeated rota- 
tions through the plume at constant axial position showed no change in the 
probe signals indicating that probe contamination was not significant. 

Hall probes . - The current distribution was measured by rotating a Hall 
effect transducer through the plume. Assuming an azimuthally symmetric current 
distribution, the self-induced magnetic field distribution is related to the 
enclosed current, J enc , b Y 


B e (r) 


^enc (r) 


2irr 


(3) 


where the assumption of symmetry can be checked a posteriori. As shown in fi- 
gures 2 and 5, the probe orientation during the traverse is such that its 
response is given by 


Bp = B r sine + Bq cos© (4) 

where B p is the field measured by the probe, B r is the radial field compo- 
nent due to the applied field and e is the angle between the probe and the 
local radius, r p . From figure 5 it is apparent that e is equal to one-half 
the angle between the probe arm and the local vertical. The contribution from 
B r is largest at the largest radius (where both 9 and B r are largest). 
Assuming current does not extend radially outside the thruster (which can again 
be checked a posteriori), the maximum radius is the thruster radius, r^, so 
that 



where R p is the length of the probe arm. For the dimensions given in fig- 
ures 2 and 5 this results in the inequality 

B r < 4B e (1) 

for the B r term to contribute less than 10 percent of the measured value. To 
establish minimum measurable currents the radial field distribution was calcu- 
lated using a thick, finite length solenoid code (ref. 18). Typical results 
are shown in figure 6 for z = 2.7 cm and three magnet current levels. For 
each magnetic field strength the minimum measurable B@ is established using 
equation (6). 

The dependence of Hall probe output on temperature is shown in figure 7. 

A 50 °C temperature increase decreased the sensor output by ~8 percent. To 
minimize the impact of probe heating during the plume traverse a 0.5-mm-thick 
mica sheet was placed on either side of the probe and the edges were sealed 
with ceramic cement. In addition, a thermocouple was located immediately 
beside the probe to monitor its temperature during thruster operation. Prelim- 
inary testing revealed that while the probe temperature did not increase sig- 
nificantly during the passage through the plume, prolonged exposure, even while 
far outside the main plume, increased its temperature by over 100 °C. For this 
reason a small water cooling line was passed up through the probe arm. High 
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thermal conductivity epoxy was used to conduct heat from the probe back to the 
arm. In this fashion the sensor temperature was kept constant without sacri- 
ficing the smail probe size. 

Fi 1 ter photography . - A CID camera was used with narrow bandpass optical 
interference filters to study global plasma species distributions. The camera 
was mounted perpendicularly to the thruster axis to image the near field plume. 
The distribution of the 4p ionized argon quantum state was examined by observ- 
ing the 4879 A line emission with a filter bandpass of 10 A. Images were 
directly digitized using a frame grabber and commercial software. At this 
stage of the research no Abel-inversion was performed. 


PRELIMINARY RESULTS 

The MPD thrusters were operated at current levels between 1000 and 1500 A 
at an argon mass flow rate of 0.16 g/s. The solenoid current was varied 
between 400 and 1500 A. Voltage-current characteristics and performance meas- 
urements are reported in reference 10. 

Electron density and temperature profiles for thrust F at several axial 
positions downstream of the exit plane are shown in figures 8(a) and (b). 

These data were obtained at a discharge current, Jt, of 1500 A, magnet current, 
Jm, of 800 A and a mass flow rate of 0.16 g/s (J 2 /m - 1.4x10^® A 2 *s/kg). The 
triple probe technique fails for electron number densities below 10^ m - 3 
because of probe sheath interactions at low densities. It is clear that while 
electron temperature exhibited strong radial gradients near the exit plane 
(reaching a peak over 2.5 eV on centerline), its axial variation was gradual. 
The electron density, however, decayed rapidly in both radial and axial direc- 
tions. Shown in figure 9 is a plot of the centerline electron density versus 
axial position and the associated least-squares curve fit. Starting from 
~lxl0 2 ® m - 3, the decay is exponential at a constant rate over the first 30 cm 
of the plume. These results are similar to those obtained by Arakawa 
(ref. 19) using a permanent magnet applied field with a 10 kW thruster 
(J 2 /m = 4.4x10'® A 2 *s/kg), though the axial decay rates observed in this work 
are lower. The different decay rates are probably due to the different shapes 
of the applied fields used in the two experiments. 

The Influence of the applied magnetic field on the density and temperature 
at z = 22 cm is shown in figures 11(a) and (b). The applied field strength 
had little impact on the electron temperature at this axial position, but the 
peak el ectron>dens i ty decreased by ~30 percent and its radial extent increased 
slightly at lower field strengths. The ratio of the gas to magnetic pressure, 
defined by 


|J = ? \ I 

B 2 /2p Q 

was evaluated using the calculated axial field strengths and assuming = N e 
and T e = T^. Results for a solenoid current of 800 A are shown in figure 10. 
3 increased exponentially with axial distance from 0.02 at z = 1 cm to 0.35 
at z = 30 cm. It is apparent from these results that the magnetic field sig- 
nificantly confined the plasma in the near field expansion. 
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The intensity decay rates for the 4879 A, Aril line for solenoid current 
levels of 1500 and 800 A are shown in figure 12. For both cases the detector 
was saturated at an intensity of 255 near the exit plane. Increasing the mag- 
netic field strength increased the intensity of the emission at a given axial 
location, indicating a corresponding increase in the population of the upper 
excited state. Least-squares curve fits reveal that the axial intensity, which 
is directly proportional to the density of the upper excited state, decayed 
exponentially for both magnetic field strengths, though there was a slight 
increase in rate for z > 3 cm. The least-squares fits were very close for 
both magnetic field strengths, with the equation 

I = I Q ex p(f^ 3 ) (z cm) (8) 

describing both curves. From this result it is apparent that the characteris- 
tic length for decay was ~2.3 cm. Using the radiative lifetime for this state 
(ref. 20) (1.3xl0 -8 s) and an exhaust velocity of ~10 4 m/s (justified by thrust 
measurements (ref. 8)), the expected decay length is 0.01 cm, which is two 
orders-of-magni tude smaller than the observed length. These estimates indicate 
there is either significant collisional excitation or Arm - Aril recombina- 
tion during the expansion, especially since the calculation does not account 
for radial plume expansion. The recombination rate was evaluated using the 
cross sections given by Drawin (ref. 21) and was found to be several orders-of- 
magnitude too small to explain the data. It is thus expected that electronic 
excitation plays a dominant role in the near field plume. 

Results of the Hall probe enclosed current measurements for thruster geom- 
etry D are shown in figures 13(a) and (b) for z = 2.7 cm for two magnet cur- 
rents at a discharge current of 1500 A and argon flow rate of 0.16 g/s. The 
maximum enclosed current occurred at a radius equal to the anode exit radius, 
indicating most of the current blown into the plume attached to the face of 
the anode. Comparison of the two figures indicates that the fractional current 
blown downstream of this axial position increases with magnet current. For 
magnet and thruster current equal to 1500 A, over 30 percent of the discharge 
current was blown downstream of 3 cm. 

These measurements have helped explain the temporal variation of the volt- 
age observed under certain operating conditions (ref. 8). For these cases the 
voltage increased monotoni cal 1y with time with no change in thruster current, 
mass flow rate, or applied field strength. As shown in figure 14, the Hall 
probe measurements clearly show that a change in current distribution was asso- 
ciated with the change in voltage with more current blown downstream when the 
thruster operates at a higher voltage. The increased path length for the cur- 
rent streamlines will increase the discharge voltage in agreement with the 
observations. While the causes of the current redistribution are still 
unknown, these measurements indicate that the voltage increase was not associ- 
ated with an electrode fall voltage. 

The results presented in this section indicate that a significant amount 
of energy transfer occurs in the plume. The fact that a significant fraction 
of the thruster current was blown downstream implies not only that there was 
ohmic heating in the plume, but also that the interaction volume between the 
applied field and the discharge current had increased. In addition, the 
applied magnetic field appears to confine the plume expansion which, by 
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decreasing axial electron density gradients, increases the collision frequen- 
cies and enhances electron-ion equilibration. In fact, the electron-ion equil- 
ibration time for the measured densities and temperatures is ~10 -5 sec, which 
for a velocity of 10^ m/s yields an equilibration length of 10 cm which is com- 
parable to the length scales of the density decay. 


SUMMARY 

An MPD thruster plume diagnostics facility, including electrostatic probes 
for electron density and temperature measurements, Hall probes for current dis- 
tribution mapping, and an imaging system to establish global plasma species 
distributions has been built and is undergoing testing. Preliminary measure- 
ments indicate that for applied field MPD thrusters operated at powers between 
30 and 60 kW the electron density distribution is strongly peaked along the 
centerline and decays exponentially with axial distance from the thruster exit 
plane. The electron temperature distribution is more uniform, dropping from 
-2 . 5 eV at the exit to ~1.5 eV at a distance 30 cm downstream. These results, 
combined with calculated applied field distributions, indicate significant 
plasma confinement and energy transfer in the plume. 
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TABLE I. - M PD THRUSTER DIMENSIONS 


Thruster 
desi gnations 

Anode channel 
diameter, 
mm 

Channel 
1 ength , a 
cm 

Anode flare 

Hal f angl e, 
deg 

Length , 
cm 

D 

26 

Jm . 2 

41 

1.1 

E 

c 26 

1 b 5.4 

10.5 

5.4 

F 

c 26 

d 5.4 

10.5 

5.4 


a Upstream end of cathode to anode exit plane. 
b Cathode diameter 1.27 cm, cathode length 2.2 cm. 
“-Inlet diameter, 

d Cathode diameter 1.27 cm, cathode length 3.8 cm. 
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FIGURE 1. - SCHEMATIC OF THRUSTER SUPPORT AND THRUST STAND. 
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PROPELLANT 

INSULATOR INJECTOR 



CHANNEL 
' LENGTH 


(A) CONFIGURATION D. 


INSULATOR 



_CHANNEL 
K “* LENGTH 

(B) CONFIGURATIONS E AND F. 



HALL PROBE 


ELECTROSTATIC PROBE 


I THRUSTER 


ROTARY MOTOR 


LINEAR ACTUATOR 


[C-89-Q6378 


FIGURE 2. - MPD THRUSTER CONFIGURATIONS (SEE TABLE 1). 


FIGURE 3. - PROBES AND MOTION SYSTEM. 
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THRUSTER 



FI6URE i|. - TRIPLE PROBE CIRCUIT FIGURE 5. - GEOMETRIC CONFIGURATION OF HALL PROBE 

SCHEMATIC. DIAGNOSTIC. 



RADIUS, CM 

FIGURE 6. - CALCULATED RADIAL MAGNETIC FIELD DISTRIBUTION 
AT THE EXIT PLANE FOR 3 SOLENOID CURRENTS, z = 2.7 CM. 


FIGURE 7. - HALL PROBE OUTPUT VOLTAGE DEPENDENCE ON TEM- 
PERATURE, B * 0.3 TESLA. 
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RADIUS, CM 

(a) ELECTRON TEMPERATURE. 




FIGURE 8. - PROFILES FOR SEVERAL AXIAL POSITIONS. 
(J M = 800 AMPS, J T = 1000 amps, M = 0.16 g/s). 


FIGURE 9. - CENTERLINE ELECTRON DENSITY VERSUS AXIAL PO- 
SITION. (J M = 800 amps, J T * 1000 amps, m = 0.16 g/s). 



FIGURE 10. - RATIO OF THERMAL TO MA6NETIC PRESSURE ALONG 
PLUME AXIS, DATA FROM FIGURE 9. 
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